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Abstract. In order to quantitatively test the ability of averaged inhomogeneous cos- 
[ mologies to correctly describe observations of the large scale properties of the Universe, 

J> ' we introduce a smoothed template metric corresponding to a constant spatial curva- 

ture model at any time, but with an evolving curvature parameter. This metric is used 
to compute quantities along an approximate effective lightcone of the averaged model 
of the Universe. As opposed to the standard Friedmann model, we parameterize this 
QQ , template metric by exact scaling properties of an averaged inhomogeneous cosmology, 

' and we also motivate this form of the metric by results on a geometrical smoothing 

■ of inhomogeneous cosmological hypersurfaces. The purpose of the paper is not to 

' demonstrate that the backreaction effect is actually responsible for the Dark Energy 

phenomenon by explicitly calculating the effect from a local model of the geometry and 
. the distribution of matter, but rather to propose a way to deal with observations in the 

I backreaction context, and to understand what kind of generic properties have to hold in 

order for a backreaction model to explain the observed features of the Universe on large 
scales. We test our hypothesis for the template metric against supernova data and the 
position of the CMB peaks, and infer the goodness-of-fit and parameter uncertainties. 
We find that averaged inhomogeneous models can reproduce the observations without 
requiring an additional Dark Energy component (though a volume acceleration is still 
needed), and that current data do not disfavour our main assumption on the effective 
lightcone structure. We also show that the experimental uncertainties on the angular 
diameter distance and the Hubble parameter from Baryon Acoustic Oscillations mea- 
surements - forseen in future surveys like the proposed EUCLID satellite project - are 
sufficiently small to distinguish between a FLRW template geometry and the template 
geometry with consistently evolving curvature. 
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1. Introduction 

On the very large scales the Universe appears to be close to a homogeneous and isotropic 
state. This is usually modeled by a locally isotropic and hence homogeneous solution 
of Einstein's equations, namely the standard Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metric. The observations of tiny temperature fluctuations of the Cosmic 
Microwave Background (CMB) radiation suggest that in the Early Universe deviations 
from this 'background cosmology' were very small, thus motivating the use of linear 
perturbation theory about the FLRW solution. The experimental confirmation of the 
predicted baryon acoustic oscillations in the CMB power spectrum [H [2], as well as 
the distribution of galaxies and clusters on the large scales [3] provide a certain body 
of evidence in favor of the 'concordance model' which results from this perturbative 
approach (see, however, [1], [5]). This standard scenario relies on the assumption that 
the FLRW cosmology correctly describes the 'background cosmology', i.e. the averaged 
inhomogeneous Universe, at all times. Even though this hypothesis may be valid in 
the Early Universe, does it continue to hold even when the Universe becomes more and 
more structured at late times? 

Answering such a question has become even more important in the light of the 
still unexplained Dark Energy phenomenon in the context of the FLRW paradigm. 
The luminosity distance measurements to type la supernova (SN-Ia) standard candles, 
when analyzed within the framework of the FLRW Universe, provide strong evidence 
for a missing component characterized by a negative pressure which, by inducing an 
accelerated phase of expansion, would be responsible for the observed dimming of far 
distant SN la (O [7]). The simplest scenario to account for these observations is a 
positive cosmological constant in Einstein's equations. This is often assumed to describe 
the energy contribution of quantum vacuum fluctuations. Nevertheless, because of the 
huge discrepancy between the particle physics expected value and the observed one, 
several alternative scenarios have been investigated. For instance phenomenological 
models such as a late time slow rolling scalar field (see reviews [HI, [9]) or the Chaplygin 
gas [To] have been proposed to describe this Dark Energy component. Alternatively, 
there have been several proposals to account for these effects through modification of 
the laws of gravitation (e.g. braneworlds [Hj, scalar-tensor gravity [12], higher-order 
gravitational theories [Hllll], AWE [T5|fTB]). 

Recently, a third alternative has been considered [T71 [T8] that aims at explaining 
Dark Energy as an effect caused by inhomogeneities. However, most of the approaches 
which include the effect of inhomogeneities still rely on the postulate that the FLRW 
solution reliably describes the effective (average) evolution of an inhomogeneous 
cosmology. For instance this is the case of inhomogeneous universe models in which 
distances are computed using perturbation theory about a FLRW background [T9t [20]. 
Other efforts abandon the FLRW model and instead restrict inhomogeneities by strong 
symmetries, employing exact solutions to Einstein's field equations like the Lemaitre- 
Tolman-Bondi (LTB) metric [2ll[22l[23l[22[25l[26l[271[28l[29l[3a 
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In this work we shall test a different approach and exploit the key-insight that 
the (large-scale) kinematics of a homogeneous-isotropic state does not necessarily 
follow the kinematics of a homogeneous-isotropic solution, especially at late epochs 
characterized by the presence of large matter inhomogeneities. Indeed, the analysis 
of backreaction effects due to inhomogeneities suggests that there is a wider class of 
(large-scale) homogeneous-(almost-) isotropic cosmological models, while smaller scales 
feature strong inhomogeneities and anisotropics that both are known to exist. In such 
a case it is natural to ask whether the emergence of Dark Energy can result from the 
breakdown of the underlying assumption associated with the FLRW cosmology. 

If no assumption on the nature of the inhomogeneities is made, i.e. if we do 
not restrict them to be small deviations from a FLRW background or obeying strong 
symmetry restrictions, we can still look at effective (average) properties of Einstein's 
equations. In the simplest case such a programme can be realized by foliating spacetime 
into fiow-orthogonal hypersurfaces, restricting the matter model to 'dust', and spatially 
averaging the scalar parts of Einstein's equations with respect to a collection of free- 
falling observers (generalized fundamental observers). Whereas such a formalism and the 
dynamical equations that govern the behavior of the averaged inhomogeneous universe 
model are well-established [34j , the explicit geometry, which lies at the basis of how we 
measure distances, is left unspecified. 

Here we suggest, as a next step, to complement the general kinematical properties 
of an averaged universe model with an explicit form of a template metric that retains 
the main properties of the standard model of cosmology, such as its isotropy and 
homogeneity on large scales, but allows for structuring on small scales. The shift in 
emphasis is from postulating a strong cosmological principle that assumes local isotropy 
about every point and hence homogeneity on all scales, to a weak cosmological principle 
that only assumes (quasi-) isotropy and homogeneity on the largest observable scales. 
In this context, we retain by assumption the usual description of the Universe at early 
times, up to decoupling. However, at late times we will have to modify this description. 

To summarize let us list our main goals and assumptions. The purpose of this 
paper is to study the influence of a geometrical effect induced by the coupling between 
backreaction and averaged spatial curvature, on top of the well-studied kinematical 
effect of backreaction on the evolution of the effective volume scale factor [M]. To realize 
this study, we propose an ansatz for the effective metric of the large-scale homogeneous 
model that is motivated by previous results on the smoothing of Riemannian metrics 
by the Ricci flow. This is used to define an effective background on which the photons 
propagate; such a background can be considered as a first refinement of the usual FLRW 
background geometry. In order to get an insight into the effects associated with our 
prescription, we consider a specific example of backreaction, namely a power law of the 
effective scale factor. Here we want to stress the fact that our aim is not to show that 
backreaction effects can be fully responsible for the Dark Energy phenomenon. Instead, 
we address the converse problem: what is necessary and what kind of generic effects are 
expected for a backreaction model to be consistent with the cosmological observations? 
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And more specifically, we are interested in understanding what specificities can allow to 
distinguish between FLRW and averaged models. 

The central concept underlying our investigation is that, although the 3-Ricci 
curvature distribution of an inhomogeneous cosmological slice can be smoothed at any 
time into a constant curvature, the dynamical evolution of the averaged curvature can 
differ from the evolution of a constant-curvature (homogeneous) model. This deviation 
has recently been quantified in the framework of perturbation theory [35l [36l [37] (see 
also [38j for an estimation in the conformal Newtonian gauge), and since we have 
arguments why a non-perturbative treatment is necessary for the effects of interest, 
we shall consider the dominant perturbative mode within a general class of scaling 
solutions to a backreaction-driven cosmology. Of course, the scaling solutions cannot 
be expected to fully represent the realistic backreaction effect throughout the whole 
history of the Universe, but it is considered here for reasons of clarity and simplicity 
to illustrate the kind of effects expected from the non-trivial geometry, in analogy to 
studies using parameterizations of the equation of state for Dark Energy. 

The paper is organized as follows. We introduce the backreaction context, the 
key-equations and free parameters of an averaged cosmological model in Section 2. 
In Section 3 we develop the ansatz for the metric, which is inspired by the study of 
the Ricci flow deformation of three-dimensional Riemannian initial data sets, and we 
use this effective metric to compute quantities along an approximate past lightcone, in 
particular the luminosity distances of cosmological objects. This template metric can 
be considered as an improvement over the FLRW ansatz that considers that photons 
follow the null geodesies of a locally homogeneous and isotropic metric. In Section 4 we 
discuss the constraints on the model parameters as inferred from SN la data and the 
multipoles of the CMB peaks and dips. In particular we calculate the cosmic distances 
in the template metric and determine the 'best-fit' models with the help of exact scaling 
solutions to the backreaction problem, which include the leading perturbative mode. 
We conclude this section with a discussion on possible tests of our main assumptions. 
Finally, we summarize the results of the paper and present an outlook in Section 5. 

2. Backreaction, effective equations and their parameters 

2.1. The backreaction context 

Let us first recall some of the main points of kinematical backreaction as a candidate for 
Dark Energy. Kinematical backreaction comprises the averaged effect of inhomogeneities 
in matter and geometry, which lead to deviations of the kinematics of the averaged model 
from that of the standard FLRW model. The existence of such deviations challenges 
the basic conjecture of the standard model, namely that the Universe is described by 
a homogeneous-isotropic solution on average. A more general approach should encode 
deviations from the standard model through kinematical backreaction terms (see, e.g. 
[39| [M] for summaries of the basic framework set out in [iQl HH |32] and discussions 
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therein). The simplest account for this effect has recently been summarized in terms 
of a morphon field, namely a scalar field that provides an effective description of the 
backreaction effects within a standard Friedmannian framework [13]. More specifically, 
the price to pay of using Friedmann's equations to describe the kinematics of a spatially 
averaged model, is the presence of an extra component. In the FLRW context such a 
component would be attributed to an exotic source (either a scalar field or a modification 
of gravity) exhibiting (phantom-) quintessence like properties [8]. 

2.2. Exact kinematics of the volume scale factor 

Let us consider a Universe filled with an irrotational perfect fluid of dust-matter with 
energy density g. We foliate space-time using the ADM procedure [H], and restrict 
our analysis to the simple case where the lapse function is constant and the shift 
vector vanishes. This choice of local coordinates is equivalent to consider a family 
of observers comoving everywhere with the fluid. We are then left with a set of spatial 
hypersurfaces of 3-metric gij{t,X^), parameterized by a universal time coordinate t, 
where (-^*)ie{i2 3} Gaussian normal coordinates on the hypersurfaces associated 
with the special choice of comoving observers. Following [ID] we define an effective scale 
factor over a compact, rest mass-preserving domain P that is contained within a spatial 
hypersurface orthogonal to the flow of the dust fluid, 

MO^(¥>V", (1) 



where Vx)^ = \T>i\ denotes the initial value of the volume. We also introduce an averaging 
operator that acts on any scalar function ^ as follows: 

■■=^ [ ^Jd'^ , (2) 
yv Jv 

with the domains volume 

Vv:= [ Jd'X , (3) 
Jv 

and J the square root of the 3-metric determinant in the spatial hypersurfaces. Then, 
by averaging the Raychaudhuri equation, we obtain a second-order differential equation 
for the volume scale factor: 

S^ + 4j,G{g)^-A = Qv, (4) 



where the extra term is the kinematical backreaction term, 

^ 2 (J/)^ - I {1)1 = l{ie- {e)^f)^ - 2 (a^>^ . (5) 



1 = 6 and // = l/2[ (^^ — Q^jQ-'- ] denote the principal scalar invariants of the expansion 



tensor (here this tensor is symmetric and defined as minus the extrinsic curvature). 
We have written Eq. in terms of kinematical scalars through the decomposition 
9* J- = ^S^j9 + a^j, where 9 is the expansion rate, aij the shear tensor, cr^ = l/2a'^^aij the 
rate of shear (summation over repeated indices is understood). 
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Averaging the Hamiltonian constraint we obtain a further equation: 

3(^)^8.G(.),-A = -<5>^. (6) 

where IZ is the 3-Ricci scalar of the spatial hypersurfaces. This equation is compatible 
with Eq. (jlj), provided that the following integrability condition is satisfied: 

\d,{Q^a%) + \d,{{n)^al) = 0. (7) 

Equation ([7]) explicitly shows the generic curvature- fluctuation coupling, i.e. the 
dynamical coupling between the averaged 3-Ricci curvature and the kinematical 
backreaction term, and it is a genuinely relativistic effect (see [M] for details). Indeed 
it is in such a fully relativistic, non-perturbative and background-free context where 
backreaction models have to be studied, if one aims to access their complete properties. 

2. 3. Effective cosmological parameters 

In addition to the volume scale factor a© we introduce the volume Hubble functional 
Hx> := dx)/ax), together with the following adimensional average characteristics: 

fi^:=g(^^)^ ; ^?:=^ ; (8) 



Thus, Eq. dHI) assumes the 'cosmic quartet relation': 

nZ + nl + nl + nl = i . (lo) 

In the following we shall also refer to the components that are not present in the 
Friedmannian context as: 

:=nl + nl . (ii) 

If = 0) then comprises the Dark Energy contribution, usually dubbed X— matter; 
here, X— matter is decomposed into its physical subcomponents. 

Finally, we introduce an effective volume deceleration parameter. 



i''^-^^ = l^Z + 2nZ-nl. (12) 



2.4- Correspondence with scalar fields: the morphon field 



It has been shown in |43]| that the backreaction effect on the kinematics can be 
mapped into an effective scalar field component <l>x), dubbed "morphon", with self- 
interaction potential t/D($x>), such that the backreaction terms comply with the 
following correspondence: 

^ ;Qv = e^l-Uv ; --i-(7^)^ = 3f/^ , (13) 
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where e = ±1 is a free parameter which allows the mapping to standard or phantom- 
like scalar fields. In this scheme, we may notice that the effective potential of the 
morphon field is directly linked to the averaged scalar curvature of space. Using this 
correspondence, Eq. (j4]), ([6]) and ([7]) read: 

3— = - A-nG ( {q)^ + qI + ?,pI)+K (14) 

K^) =8vrG'((f?)p + f?f)+A (15) 

+ ^Hv^v = - ^T^-Ui^v) , (16) 

where = 6^$!, + f/p and = ei<l>|, — f/p are the energy density and pressure of the 
morphon field, respectively. This system shows that, at the kinematical level, i.e. as far 
as the evolution of the effective scale factor is concerned, the backreaction effect and a 
scalar field component cannot be discriminated. Nevertheless, the backreaction effects 
are partly of geometrical origin, and in the following we shall see that such a feature 
is the key point that, through the analysis of observations, can discriminate between 
backreaction and a standard quintessence component. 

Overall, our model can be summarized by the following argument: the proposed 
model is kinematically equivalent (up to the involved dependence on the averaging 
scale) with the standard model supplemented by a quintessence field. The difference 
is, besides the obvious physical justification of the quintessence field, dynamical: while 
the standard model has only physical sense as a template for the average distribution 
(since the Universe is not homogeneous), we have to make sure that the exact integral 
properties are satisfied as well. Hence the alternative model that we propose here can 
be considered as an attempt to 'repair' the standard FLRW scenario, since this latter 
is not compatible with the underlying metrical properties of the low redshift matter 
distribution, while the former is. Nevertheless, our model maintains some of the naive 
assumptions of the FLRW case, notably the fact that that hght propagates in a constant 
curvature space, while in reality it travels, most of the time in 4-Ricci fiat regions due 
to the volume-dominance of voids at moderate redshifts [15] . 

3. Effective geometry 

3.1. Motivation for the effective metric: construction of distance measurements 

Our approach is the following: 

i) Space-time is foliated using a 3+1 ADM procedure [44j- Thus the metric reads: 

(^)g = -dt^ + gijit, x'')dX' ® dX^ , (17) 

where t is the cosmic time, latin indices are in {1, 2, 3} and gij{t, X^) is the metric 
of the spatial hypersurface at time t (first fundamental form), endowed with the 
local comoving coordinates X'^. 
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ii) An effective averaged model is then built by the averaging procedure described in 
|40j . [H] and |12]. This results in the Eqs. (jlj) and ([6]) which contains the so- 
called backreaction terms stemming from the averaged fluctuations in the local 
expansion and shear scalars. These equations admit solutions for the effective scale 
factor avit), which describes the kinematics of the averaged model, provided that 
the integrability condition Eq. ([7]) is satisfied. We shall close this system by the 
ansatz of exact scaling solutions |13] with the scaling index as a free parameter. A 
particular scaling index reproduces the result of second-order perturbation theory 
in the synchronous gauge [35] . 

iii) In addition to accounting for the kinematical properties of backreaction models a 
definition of distance measurements is still required for a complete description. 
So far, the structure of the lightcone has only been explored by perturbative 
calculations, e.g. [HI [191 120], but our modification of the standard model 
focuses on the effective cosmological evolution and the determination of effective 
distances within a backreaction model that, without any assumptions on the 
kinematical properties, couples the evolution of fluctuations to that of the averaged 
scalar curvature. This is the reason why our modification is inspired by recent 
results that deal with a rescaling of inhomogeneous three-geometries as a possible 
solution to the so-called fitting problem [171 1^ . Ricci flow renormalization of the 
average characteristics on a bumpy geometry (see [HI HI] and references therein) 
monotonically decreases intrinsic curvature inhomogeneities and would produce a 
constant-curvature slice, but only at a given instant in tim^. This result then 
provides us with a prescription to construct distance measurements along the 
lightcone. As an effective model for the geometry of our Universe, we choose to 
replace, at a given scale corresponding to the domain "D, the true 3+1 foliation 
given by the metric, Eq. ( |T7I) . with another 3+1 foliation such that at each time 
the spatial hypersurface is modeled by a Riemannian 3-space of spatially constant 
curvature, which we denote as K,v{t*)- Such a curvature term can obviously be a 
function of time t*, since there is no reason why an homogeneization procedure 
such as the Ricci flow, that is applied separately at any time, would result in 
the same constant curvature space. In contrast, the coupling between spatial 
fluctuations in the second fundamental form and the smoothed curvature implies 
that, if the fluctuations have changed from one instant in time to another, the 
resulting smoothed curvatures at these two instants may be different. In such 
a case a template model of space-time consists of a family of three-dimensional 
hypersurfaces of spatially-constant curvature, with each hypersurface being 
characterized by two scalars: ax)(t) for the kinematical properties and K,v{t) for 
the curvature. Therefore, the model is specified by (T,t,aT>(t), KD(t))^^^. Since 

I In general, such a flow has singularities, if the 3-Ricci tensor is non-positive, and a constant -curvature 
model is reached only after subsequent steps of surgery of the manifold. Employing surgery, Perelman 
[50l [51] has added the final proof of the possibility of smoothing a generic 3-manifold into pieces of 
constant-curvature manifolds with the help of the Ricci flow. 
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the kinematics of an averaged model is governed, in general, by the volume scale 
factor a-D and not by the Friedmannian scale factor a, introducing the former in 
the effective description of the lightcone is a justifiable assumption (see [52] ji- 

We can now compute the geometrical properties of the effective model. Let us 
choose on each hypersurface the standard spherical coordinates (r, 6',0). 

At a given time the distance between two points on a hypersurface Ej, separated 
by a coordinate distance r is given by: 

r dr' 



a/1 - KD(t*)r'2 

Consider a second hypersurface St.+dt infinitesimally close to Sf, . We can write 
the distance between the two points over this new hypersurface as 

r dr' 

Kt* + dt) = av{t, + dt) / . (19) 

Jo V 1 ~ + dt)r'^ 

By Taylor expanding Eq. f|T9l) around t^,, keeping only the first order terms in dt and 
using f[T5]) . one obtains: 

§(t) = H^mt) + av{t)'^-^{t) f '''^"^''^ . (20) 
dt dt Jq 1 — Kx)(t)r^ 

The first term on the right-hand-side is the standard Hubble law, and the second 
term expresses a deviation from this Hubble law induced by the particular geometrical 
structure of the effective model. Let us consider now the motion of photons. Since their 
velocity is constant and equal to c along their geodesies, one has the relation: 

(21) 

Hence, using Eq. ( ITHl) . one finds: 

diT 

cdt = dlit) = ap = , (22) 

a/1 - Kv{t)r^ 

or equivalently, the following differential equation for the coordinate distance travelled 
by a photon: 

dr c 



dt ax> (t) 



1 - Kv{ty . (23) 



§ Note here, that the values of these two scale factors are not expected to differ strongly on large scales. 
On the other hand, their derivatives can differ significantly and this is accounted for by allowing for an 
explicit and consistent change of the scalar curvature. As we assume that the standard FLRW model 
describes the Early Universe, including the recombination epoch, backreaction effects are assumed to 
be suppressed and consequently both scale factors coincide at early times. 
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3.2. The template metric 

The ordinary differential equations ( !20l) and ( !23|) suggest to introduce an effective metric 
to describe the geometry in the 4-dimensional template averaged model, so that the 
photons move on null-geodesics of this metric. These equations are exactly the same as 
those that would be inferred from the following space-time metric: 

= -de + L]j^al-ff^dX'®dX^ , (24) 

where a^ioLHo = ^/Hvo (c = 1) is the present size of the horizon introduced so that the 
coordinate distance is adimensional, ax>(t) is an adimensional volume scale factor on a 
rest mass-preserving compact domain V, endowed with a domain-dependent effective 
three-metric that assumes the form: 

where dfl"^ = r'^^dO'^ + s\v?{6)d(j)'^) is the solid angle element. This form has been 
suggested, using different considerations, by Paranjape and Singh [52] . 

At any given time, the template three-metric is identical to the metric of the spatial 
part of a FLRW space-time, but with a scalar curvature that smoothly varies from time 
to time. 

However, for consistency with the general kinematical properties discussed above, 
n-D cannot be arbitrary, and it must be related to the true averaged scalar curvature 
(7?.) 25. We propose, in analogy with a FLRW metric, to choose such that: 

alit) 

This choice guarantees that K-p does not have any dimension, and suitably inherits the 
sign of {TV) J,- We shall discuss the possibility to test this prescription at the end of the 
paper. The FLRW case can be recovered by formally posing K©(to)| (^)x>o I ~ ^©0/65 
with /tx)(to) '■= 1; where fcx^o is identified with the constant Friedmannian curvature on 
a chosen large domain Vq. 

It is worth emphasizing that this template metric is not required to be a dust 
solution of Einstein's equations [53], [51] (the effective fluid of an averaged dust model 
also features a geometrical pressure [M])- In fact Einstein's field equations are satisfied 
locally for a general, unspecified space-time metric, whereas the template metric defined 
by (12^ and (!25l) is only a prescription to take into account the non-trivial behavior of 
the three-curvature in the averaged effective model of the Universe, and to compute 
quantities along the approximate smoothed lightcone associated with the travel of light 
in a clumpy Universe (for another earlier attempt see [55])- Recently, Kasai [56] has 
investigated the goodness of fit to supernovae data of Friedmannian models without 
cosmological constant, and with different curvaime parameters. While a single standard 
model without cosmological constant cannot account for the supernova data, two such 
models - if applied to low- and high-redshift data separately - can reproduce the 
observations [56]. Kasai also provides in his paper formulae from a post-Newtonian 



(TZ)^ = ^° . (26) 
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description of inhomogeneities to explain dynamically the link between these two models. 
Our (background-free) model differs in this particular aspect, since we work with general 
averaged cosmologies, e.g. we do not assume the scale factor to be that of a FLRW 
solution. Also several other works have concentrated on the role of curvature (however, 
within the class of constant-curvature models), e.g. [57], [58] . 

We shall show below that supernovae data by themselves indicate that the 
assumption fl26l) is indeed reasonable. From now on, we consider this effective metric to 
compute cosmic distances. 



3.3. Computation of observables 

The computation of distances is a very important issue that must complement the 
kinematics when comparing the predictions of a given averaged model with observations. 
In an inhomogeneous Universe, this is a difficult problem that can be overcome either 
by considering exact solutions of Einstein's equations with particular symmetries such 
as the LTB metric, or by explicitly computing the photon path through a perturbed 
FLRW background [IHl [20]. 

The above prescription considerably facilitates the computation of effective 
distances. Firstly, let us introduce formally an effective redshift zx>. The redshift is 
defined as: 

where the O and the S stand for the evaluation of the quantities, respectively, at the 
observer and at the source. In this expression, gab is the effective metric (!24l) . the 
4- velocity of the matter content {u°'Ua = —1), supposed to be comoving hereafter, and 
k"- the wave vector of a light ray travelling from the source S towards the observer 
O {k^ka = 0). Then, by normalizing this wave vector such that {k^Uajo = —1 and 
introducing the scaled vector k"" = a'^k"', one obtains the equation: 

l + zv= iai'k')s , (28) 

with k^ obeying the null geodesies equation k^'W ak^ = which, since in our case we are 
only interested in leads tc[3: 

1 dk° r'^iav) dKv{av) 



kpdav 2(1 - Kx,(ax))r2(a-p)) da-D 
As usual, the dimensionless coordinate distance can be derived from the equation 
of radial null geodesies, as was done for Eq. ( |23l) : 

Solving Eq. (I5U]) provides the coordinate distance f(ax>) that can then be plugged 
into Eq. fl29l) to find the relation between the redshift and the scale factor. With both 

II See the note at the end of the paper. 
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r(ax)) and z{ax>), one can determine the angular diameter {dA{zx>)) and the luminosity 
(c^l(^I))) distances of sources respectively, which read: 

dA{zv) = -T^av{zv)f{zv) , (31) 
dL{zv) = (1 + zvfdA{zv) , (32) 

the last relation holding because of the reciprocity relation between distances. Having 
computed these functions, it is then possible to compare the backreaction model 
predictions with the data. 

The observations of SN standard candles provide direct measurements of the 
luminosity distance as function of the redshift, thus they can be used to constrain 
the effective parameters (previously introduced) without the need for any further 
assumptions on the model. On the other hand the analysis of the measurements of 
the angular distribution of the CMB temperature fluctuations requires the solution 
of the Boltzmann equation together with the Einstein equations. So far this has 
been derived only in the context of the linear perturbation theory about the FLRW 
solution. It is beyond the scope of this paper to perform this calculation in the 
backreaction context. Nevertheless, as mentioned before, it is reasonable to assume that 
the Early Universe (until the recombination epoch) is well-described by a perturbed 
FLRW modei^. Therefore, with this assumption, we can safely limit our analysis to 
only using the cosmic distance information encoded in the location of the acoustic 
oscillations in the CMB power spectrum (see [591 ED] , and for a recent detailed analysis 
|61]). In fact the multipoles of the CMB peaks and dips can be determined from the 
measured anisotropy power spectrum without the need of making any assumption on 
the underlying cosmology of the late-time Universe, provided this cosmology is well 
approximated by a FLRW Universe in the context of General Relativity before and up 
to the time of recombination. These multipoles can then be compared with the model 
prediction given by 

Im = la{m - (pm) , (33) 

where m eN* for peaks and m G {3/2, 5/2, ...} for troughs, and 

L = TT— , (34) 

where Oxj^rec is the scale factor at the recombination, f(ax),rec) is the comoving distance 
between recombination and an observer located on Earth today (obtained by integrating 
Eq. (13T|) ). and rs(ax),rec) the sound horizon size at recombination given by multiplying the 
sound speed c^/c by the result of the integration of Eq. (1311) ) between ax>,rec and ajy = 0. 

% In fact we have no other choice until effective properties of an inhomogeneous inflationary model 
are analyzed in detail, which would help to decide whether we can identify a homogeneous-isotropic 
state at the exit epoch with a homogeneous-isotropic solution, or whether we have to deal with a 
backreaction remnant also in the initial data. We have to admit that our arguments of a coupling of 
fluctuations with geometrical properties would also apply in the inflationary context. 
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The term 0^ parameterizes the effect of dragging gravitational forces occurring before 
recombination and which displace the maxima and minima of the acoustic oscillations 
with respect to the sound horizon scale [60] . The amplitude of these shifts only depends 
on the matter and baryon densities, on the redshift of recombination, on the scalar 
spectral index for primordial fluctuations and on the density of Dark Energy at the 
time of last scattering. Parameterization formulae for 0^ are given in [601 E2] . They 
depend on the following set of parameters: the matter density Vl^ , the baryon densitjo 
Q^°h?, the spectral index n^, and the Dark Energy density at the last scattering surface. 
Since we have supposed that the Universe is almost Friedmannian at the time of last 
scattering and since in our model the density of Dark Energy is directly related to 
the deviation from a Friedmannian universe model, we shall neglect the dependence on 
Early Dark Energy. We also impose that the spectral index is 1, corresponding to 
a scale invariant primordial power spectrum. Then, the evaluation of ipm only depends 
on pre-recombination quantities that we include in our likelihood analysis. We wish 
to remark that the only quantity that is specific to our model is /a, which is explicitly 
computed for the template metric fl2^ by integrating Eq. fl30|) . la is also used in the 
alternative approach by [63], but their second quantity, i?, cannot be used safely since 
it depends on the behaviour of perturbations [61]. 

4. Constraints from supernovae and CMB observations 

We perform a likelihood analysis to infer constraints on the backreaction models using 
the SN la data from the SNLS collaboration [65| and the position of the CMB peaks 
and dips from WMAP3-yr data measured in [66]. In particular we focus our analysis 
on a simple class of backreaction solutions, the so called 'scaling solutions', that we are 
going to describe below. 

4-1. Exact scaling solutions 

Scaling solutions have been extensively studied in [13] with A = 0. Following the 
notations of [13], we look for solutions to the backreaction problem in the form: 

Qv = Qv^l ; {'^)v = (^)a «?> ' (35) 

where n and p are real numbers. In [43j, two types of solutions were found. The 
first type, with n = —2 and p = —6, is not very interesting for our purpose, since 
at late time it corresponds to a quasi-Friedmannian model in which backreaction is 
completely negligible. Kinematical backreaction and averaged scalar curvature decouple 
in this degenerate situation. This situation corresponds to a 'quasi-Newtonian' evolution 
of fiuctuations on a non-dynamical background geometry. In contrast, the second 
type of solutions corresponds to a direct coupling between kinematical backreaction 
and averaged scalar curvature, for which we find n = p, and Qt> = r {TZ)jy with 

+ This quantity is defined exactly as in and h := H-p^/lOO as usual. 
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r = — (n + 2)/(?2 + 6). Henceforth, we shall consider backreaction and averaged scalar 
curvature of the following form, just retaining the generic class of scaling solutions: 

(7^)^ = (7^)^^ al (36) 

Qv = - (^) A «S> • (37) 



n 



Defining 



one has: 



- ~wTmi ' ^^^^ 

i/|,(a^,) = Hl^ {nZ^a^^ + r^^^a^) (39) 
t + 6)r]^°a^" 

T(^rr6j?]|^ 



. _ . (40) 



Then, Eq ( l30l) becomes: 



- / ' - ; r(0) = . (41) 



The leading order of the perturbative estimation of the backreaction effect (i.e 
the leading order at late time) is included in this class of solutions, and corresponds 
to the case n = — 1 [3S]. This can be considered as the minimal estimation of the 
full backreaction effect, since it accounts for the backreaction in a quasi-Friedmannian 
Universe, thus very close to a locally homogeneous and isotropic configuration. For 
simplicity we assume the full backreaction effect to be described by a power law 
solution, although a realistic treatment of the backreaction effect will not correspond to 
a fixed scaling, since we expect a qualitative change at the epoch of nonlinear structure 
formation. 

In the Dark Energy context, the kinematics described by Eq. ( 139|) corresponds to 
introducing a component with constant equation of state 

u;^ = -i(n + 3). (42) 

However, we must not forget that the geometry is different due to the time-dependence 
of the curvature. The leading order perturbative solution corresponds to w = —2/3, 
while conversely a cosmological constant (that is a particular morphon and not explicitly 
introduced) would require n = 0. From the scaling solutions previously described, we 
can reconstruct the morphon field; its potential during the matter dominated era is 
given by [13]: 

" 24^^G V ) 

2(l + r) 




((1 + rh^^) (^^) A sinhl^ ^v^$^ , (43) 
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where {Q)vi is the averaged restmass density of dust matter at a given initial time, and 
7^j^ = fi^'/fi^' is the ratio of energy density in averaged curvature and dust at the same 
initial time. This potential is well-known in the quintessence context as it corresponds 
to quintessence fields with a constant equation of state, or in other words, to a constant 
fraction of kinetic and potential energies (63, EH 169] . 

This equivalence shows that, at the level of the dynamical equations for the scale 
factor ax>i the effect of backreaction cannot be distinguished from the effect of a 
minimally coupled scalar field in a FLRW context. Nevertheless, the geometrical effects 
induced by backreaction along the lightcone, and introduced in Subsections 3.1 and 
3.2, are a particular feature of a backreaction model and can be a way to distinguish 
between a FLRW model and a backreaction model. To illustrate this fact. Fig. [1] shows 
the evolution of the coordinate distance cr/H-p^ in two models with the same parameters 
n and one corresponding to a FLRW Universe, and the other one to a backreaction 
model with a compatible geometry of the lightcone described by Eq. (125|1 and Eq. (126|) . 
It is clear that distance measures like standard rulers or standard candles can distinguish 
between the two models. 

4-2. Constraints 

In Fig. [2] we plot the joint constraints with la and 2a likelihood contours obtained from 
supernovae and CMB peak locations for this class of models (filled contours); we also 
plot the case of a spatially fiat FLRW model in presence of quintessence with equation 
of state given by Eq. fj42|) (solid lines). 

It can be seen that the averaged geometry slightly pushes the constraints towards 
higher values for compared to standard Dark Energy models with a FLRW 
geometry. On the other hand, although a high level of backreaction is needed and 
acceleration of the effective scale factor still occurs, a lower amount of backreaction 
is required when accounting for the effective geometry. This is particularly evident 
regarding the absolute best-fits. The leading perturbative mode (n = —1) as calculated 
in [35] and [37], is marginally at Icr for values of ~ 0.3, whereas a Dark Energy 
model in FLRW with n = —1 is compatible with the data at la for ~ 0.1. 
However, as expected, the geometrical effect is not sufficient to reconcile the observations 
with a purely perturbative estimate of the backreaction term on horizon scales |37j ; 
nevertheless, as a first approximation of a compatible structure of the lightcone, the 
effect acts in the right direction. In Fig. [2], we have represented the absolute best- 
fit for each model: the diamond is for the averaged effective model, corresponding to 

= 0.38, n = 0.12, Ht,„ = 78.54 km/s/Mpc and nf"h'^ = 0.0255, and the disk is 
for the fiat FLRW model, with Q'g" = 0.18, n = -0.5, Hjy^ = 82.0 km/s/Mpc and 
Q^°h? = 0.0263. The difference in the between the two models is 3.9. The best-fits 
marginalized over all other parameters (i.e. ifx>o ^b""^^) given by = 0.397 
and n = 0.5 for the averaged model and by = 0.26 and n = 0.24 for the fiat FLRW 
model. The large difference between the absolute and the marginalized best-fits for the 
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Figure 1. Evolutions of the coordinate distance cr/Hx>o (panel a), the redshift (panel 
b) and the ratio between the FLRW redshift to the redshift, l/{ax>{l + z-p)) (panel c), 
in the averaged model as functions of the effective scale factor, in an averaged effective 
model with n = —1, 17^° = 0.3 and Hx)o = 70 km/s/Mpc (dashed line). In Panels a 
and b, the solid lines represent a FLRW model with the same set of parameters. Panel 
d represents l/(aD(l + z-p)) for the best-fit averaged model described in subsection 
4.2. 



averaged model parameters is caused by the non-Gaussian structure of the hkehhood 
function. Fig.[3]shows the evolution of the adimensional density parameters (right panel) 
and K-xy (left panel) respectively for the best-fit backreaction model, corresponding to 
the diamond in Fig. [2] for which n = 0.12 and = 0.38. The ratio between the 
effective redshift in this best-fit and the standard FLRW redshift, l/(ax)(l + zx>)), can 
be found on the panel d of Fig. [TJ The difference in the early epoch is of order 25 %. 
From the likelihood plot (Fig. [2]), we can see that the likelihoods for the two models 
are displaced along the direction orthogonal to the degeneracy line. For fixed values of 
n and fl^ the coordinate distance of the backreaction model is systematically smaller 
than for the standard Dark Energy model. Hence, in order to fit the data, standard 
Dark Energy models require more negative values of the equation of state (i.e. larger 
values of n) and smaller values of compared to the backreaction models. 
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Figure 2. Supernovae and CMB constraints in the (0^°,n) plane for tlie averaged 
effective model with zero Friedmannian curvature (filled ellipses) and for a standard 
flat FLRW model with a quintessence field with constant equation of state w = 
— {n + 3)/3 (black ellipses). The disk and diamond represent the absolute best-fit 
models respectively for the standard FLRW model and the averaged effective model. 

4-3. Testing the curvature prescription / f^) 

In the process of constructing the averaged effective model, we assumed Eq. ( |26|) . namely 
that the geometrical instantaneous spatially-constant curvature ^©(t) is related to the 
actual averaged scalar curvature {71)^. As previously emphasized, this prescription 
is reasonable from the physical perspective, nonetheless it is worth to test is validity 
directly against the data. Indeed, the presence of this non-trivial curvature effect, 
that makes the effective template metric compatible with the averaged scalar curvature 
(7^)p, is purely the result of inhomogeneities. Whereas from a kinematical point of view, 
acceleration driven by standard quintessence or by backreaction are indistinguishable, 
a non-trivial curvature effect favoured by the data would be an unambiguous signature 
of the presence of backreaction. Supposing that the scalings are: 

Kv{av) oc and {TZ)^ oc a^+^ , (44) 
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Figure 3. Left panel: evolution of n-p for the best-fit averaged effective model 
(diamond in Fig. [2]) with n = 0.12, ^1^° = 0.38, and Ha 78.54 km/s/Mpc. Right 
panel: evolution of il^ (dashed line) and il^ (solid line) for the same model as in the 
left panel. 



where n and m are arbitrary real numbers, we shall say that the assumption ( !26ll is 
valid if and only if n = m. Remember, that a pure scaling ansatz is not what we expect 
in a realistic evolution of backreaction, so we would be surprised if the data favour this 
relation exactly. 

Figure m shows the marginalized likelihood for the variable n — m, the vertical line 
represents the case n = m. As we can see, assumption ( |26ll lies within la of the best-fit; 
it is consistent with current data, but the constraints are not very strong and currently 
depend on prior limits for m and n (here we choose both of them to lie in between -4 
and 4). It will be interesting to perform this test when more accurate data become 
available, as well as accounting for better closure conditions of the averaged Einstein 
equations rather than a simple scaling solution, hence describing the time-evolution of 
backreaction more realistically. 

4.4- Measuring the evolution of curvature directly 

Another promising way to test Assumption ( l26l) relies on a quantity introduced recently 
in [To] to test the Copernican principle. In a constant curvature FLRW universe, the 
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Figure 4. Probability distribution function for the difference of scalings n — m, 
marginalized over all other parameters. The vertical line indicates the case n ~ m 
corresponding to Assumption ([26]) . 



coordinate distance D = r / Hq is given by: 

1 



HnD z) 



: sm 



k,0 



Hiu) 



-du 



(45) 



Taking the derivative of r{z) with respect to z (denoted by a prime in the following) 



turns the sine into a cosine, and it is easy to see that {H{z)r\z)Y 
for any redshift z. The quantity 

{H{z)D'{z)f - 1 



-QkfiiHoriz))' 



(46) 



is therefore constant and equal to Qk,o in a FLRW universe. Replacing r{z) by f{zv), 
and considering the solution to equation ( l30l) integrated between and zj), we find that 
Qf{zTi) is not constant, but equal to —kv as can be seen from Eq. (130|) . 

Evaluating the derivative of flk{z), after some algebra we derive the function C{z) 
defined by equation (3) in [70]. This function identically vanishes in a FLRW model, 
whereas for our lightcone prescription it reads: 

Hv(zv)f^(zv)i^'vizv) 



(47) 



It is evident that, if the effective geometry of the Universe reduces to FLRW, then kd is 
constant, and Cx> vanishes for all redshifts, as stated in [7D]. In Fig. Owe show Cx> for 
the best-fit scaling solution corresponding to the diamond in Fig. [21 that is for n = 0.12 
and = 0.38. There are some features common to all the averaged effective models 
based on scaling-backreaction that fit the SNla and CMB data: they all have a Cx> 
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close to zero at the recombination redshift (corresponding to our assumption that the 
Universe is well described by FLRW at that time), and a decrease in the late-time 
Universe, with a minimum at a redshift around 3 to 4. 




1 10 100 1000 

l + z„ 



Figure 5. Evolution of Ct> as a function of 1 + 2-d for an averaged effective 
model with scaling backreaction whose parameters are: n — 0.12, = 0.38 and 
i/po = 78.54 km/s/Mpc (absolute best-fit model represented by a diamond in Fig. [2]). 

The form of C-p or in the averaged effective models (H7|) . and its evolution in 
the best-fit presented in Fig. O can be considered as predictions of the averaged effective 
model, that could be tested in order to discriminate it from a standard FLRW Dark 
Energy model. In particular, all models explaining Dark Energy through effects from 
structure formation will show a strong departure from FLRW at late times, z-p ~ 10. If 
this feature does appear in the observational data, it will be an unambiguous sign that 
the FLRW approximation is not valid in the late-time Universe, and if the behavior 
of Cx! follows a curve similar to the one presented in Fig. [5], it will indicate that the 
template metric employed in this paper in conjunction with the backreaction effect are, 
despite its obvious shortcomings, a satisfactory description of the late-time Universe. 

Measuring the function Cv{zv) or even VL^{zv) remains difficult because of the 
necessity of evaluating the redshift derivative of the coordinate distance directly from 
the data (see the upper panel of Fig. [6]). However, we may notice that, not surprisingly, 
all necessary information is contained in Hx>{zxi) and r(zx))- In fact, in order to determine 
the curvature, it is necessary to measure observables which depend on the geometry, such 
as the luminosity or angular diameter distance, both derived from r{zD), and another 
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observable that defines the kinematics, i.e. Hx>{zx)). Therefore, to test the geometry, 
it is sufficient to combine measurements of r{zxi) and Hxi{zt,)- If the Universe evolves 
with a constant curvature, the predictions of the averaged model for r^z-v) and H-p^zv) 
will not fit the data simultaneously; vice versa, if the curvature is not constant, then it 
will be the FLRW model that badly reproduces the observations. Unfortunately there 
are currently no strong experimental limits on the evolution of the Hubble parameter. 
To get an idea of the size of the effect, we will assume that we can use Baryonic 
Acoustic Oscillations (BAO) to measure both Hd{zv) (from the radial component) 
and Da{zti) = f{zv)/{^ + Zt>) (from the transversal component), although there may 
be some complications (see [70j for more discussions). As we shall see shortly, the effect 
is small and difficult to measure, so that we consider a future survey extending over a 
large redshift range, similar to the one expected for the EUCLID satellite projecO This 
survey would be able to measure the BAO in 9 bins in the range of zv = 0.4 to 2.0 with 
an expected accuracy of about 1% in each bin. As shown in Fig. [6l a flat FLRW with 
a standard quintessence having the same parameters as the averaged cosmology model 
can be ruled out by EUCLID: even if the two models share the same kinematics, i.e. 
the same Hubble parameter, they lead to different results for the distance measurement, 
hence they could be discriminated. 

The central and lower panels of Fig. [6] show respectively the behavior of the distance 
and Hubble parameters in the best-fit averaged model (solid line), the FLRW model 
with the same parameters (dot-dashed line), and a FLRW model with a cosmological 
constant and the following parameters (dotted line): flrn,o = 0.277, ^a^o = 0.735 and 
Ho = 73 km/s/Mpc. The latter model was chosen to reproduce the angular diameter 
distances of our best-fit averaged model, including the location of the CMB peak. Even 
if the averaged model correctly describes the Universe, there are FLRW models that 
can reproduce the distance data, but then cannot account for the Hubble rate, and vice 
versa. This confirms that EUCLID measurements together with the CMB peak position 
will be able to discriminate between FLRW models and averaged cosmologies. For a 
best-fit FLRW model as the ACDM scenario presented here, this will require the use of 
a combined measurement of the distance-redshift relation and of the Hubble parameter 
at different redshift, whereas for a FLRW model with the same set of parameters as 
the best-fit averaged model, a single measurement is sufficient. This latter point is a 
consequence of the non-trivial relation between redshift and scale factor in the averaged 
cosmology; indeed, even though the two models share the same Hubble parameter as 
a function of the scale factor (as a result of their sharing the same parameters), the 
relation between redshift and scale factor is different in the two models, then resulting 
in different Hubble rates at a given redshift. 

* EUCLID is composed of the former proposed DUNE [TT] and SPACE 75] missions. 
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Figure 6. Upper panel: Evolution of Vlk{z-D) as a function of redshift for the absolute 
best-fit averaged model represented by the diamond in Fig. O The error bars are 
obtained using the EUCLID satellite project measure uncertainties. One can see that 
all positively curved FLRW models (ilfc,o < 0) and only highly negatively curved FLRW 
models {^kfi > 0.5) can be excluded by the estimation of rifc(z-p). Central panel: 
Evolution of the coordinate distance for the best-fit averaged model (solid line), for a 
ACDM model with ri,„,o = 0.277, I^a = 0.735 and = 73 km/s/Mpc (dashed line), 
and for the FLRW model with the same parameters as the best-fit averaged model 
(dashed-dotted line). The error bars are still obtained using the EUCLID expectations. 
Lower panel: Evolution of the Hubble parameter H/Hq for the best-fit averaged model 
(solid line) , the FLRW model with the same parameters as the averaged best-fit model 
(dashed-dotted line), and for the same ACDM model as in the central panel (dashed 
line). The error bars are still obtained using the EUCLID expectations. 
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5. Conclusion and Outlook 

In this paper we have addressed the problem of comparing averaged inhomogeneous 
cosmologies with observational data by proposing to fit the observations with the help 
of an improved template metric, whose form is compatible with the kinematical integral 
properties of a general averaged model. This template metric has been motivated by the 
fact that the averaged curvature of space-like hypersurfaces of the space-time foliation 
is not expected to be constant in time. Indeed, the cosmological principle only requires 
the spacelike quantities to be averaged, but does not impose anything on the evolution 
of these quantities. In other words, the FLRW universe models, which obey a strict 
cosmological principle, are a very particular subclass of models respecting a weaker 
cosmological principle presented in this paper. We consider the modified template metric 
as a first approximation tool for interpreting observations in a Universe that appears 
homogeneous on large scales, but in which the backreaction effect cannot be neglected. 
That means that the proper effective lightcone along which the cosmological observations 
are made cannot be simply approximated by a FLRW lightcone. It is indeed important 
to notice that this template metric has only been introduced to compute quantities on 
the lightcone. 

Thanks to this prescription for the lightcone, we then deduced constraints on the 
particular class of scaling backreaction, using the supernovse luminosity /redshift 
distribution, and the positions of the peaks in the CMB spectrum. We found that the 
non-trivial geometry of the lightcone induces a slight change in the constraints, with 
respect to the same models in a FLRW geometry, leading to models compatible with 
the data for higher values of that is to say with a smaller amount of backreaction. 
This is particularly true for the leading perturbative mode (n = —1). 

One should note that the model presented in this paper still needs an acceleration 
of the effective volume scale factor to reproduce the data. Recent results based on the 
the Lemaitre-Tolman-Bondi (LTB) solution (see [32] and references therein) show that, 
on the contrary, one can fit the data without an acceleration of the volume scale factor. 
The two results do not necessarily disagree. Indeed, in the LTB model, one fits the data 
with an inhomogeneous metric having a functional degree of freedom {th{r) in [32]), and 
then, one averages the best-fit model to find that there is no acceleration of the effective 
volume scale factor. In this paper, we have first introduced an effective homogeneous 
model that we fitted to the data, and instead of having a functional dependence in 
the effective metric, we have made the assumption that backreaction features a scaling 
behavior with the scale factor. The comparison between the two approaches could 
be done once more realistic inhomogeneous metrics are at hand, and a more realistic 
behavior of the backreaction effect could be implemented. Volume acceleration might 
be a consequence of over-evolving the backreaction at early times by a strict scaling 
ansatz. 

One important feature of our results is that the cosmic history and the distances are 
strongly affected by the introduction of a non-FLRW geometry for the past lightcone. 
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even if the constraints are only slightly different in the (f2^,n) plane. 

Finally, we discussed two ways of testing the assumption fl26|) that links the 
geometry of the lightcone with the kinematical properties of the model. On the one 
hand, we found that the supernovae data are fully compatible with the assumption 
that the averaged curvature felt by photons along the past null-cone is linked with 
the averaged Ricci scalar of space-like hypersurfaces according to fl26l) . Unfortunately, 
the current available data are not sufficiently precise to unambiguously show a preferred 
selection of this hypothesis; this analysis should be done again with future data providing 
more statistics. On the other hand, we have calculated the explicit form of a function 
Cv^zo), previously introduced in [TOj to measure possible violations of the Copernican 
principle. This form can be considered as a prediction of the particular models studied 
in this paper, and it will be a crucial test demonstrating a quantitatively significant 
difference to the standard FLRW paradigm. We have shown that the future EUCLID 
satellite project might be able to distinguish between a FLRW template geometry and 
the template geometry with evolving curvature in an averaged model by using joint 
measurements of a geometrical quantity {dA^z)) and a kinematical property {H{z)). 

The effective metric that we employed has been motivated by physical and 
mathematical arguments, but it cannot be considered as a fully satisfactory description 
of the lightcone. In particular, since weak lensing involves a series of local effects of 
the metric on lightrays rather than just an integrated effect, implementing constraints 
from weak lensing surveys require a more refined study of the lightcone structure 
and of its link with spatial averaging. This link can be achieved by implementing 
the averaging formalism directly on the lightcone; this is the subject of work in 
progress. Another improvement will come from a closure condition that is better than 
a simple scaling solution, and that will encode more precisely the time evolution of 
backreaction. Such a closure condition can be looked for in numerical simulations, 
analytical approximations (see [73l[711[75] for particular approximations, and [76 11771178] 
for an interesting perspective), or observations (see [79] and references therein for 
remarks on the difficulties of this last approach). Moreover, the complete study of 
other observables like the full CMB spectrum is still unavailable and will be crucial 
for the construction and the test of a 'concordance model' for averaged inhomogeneous 
cosmologies. 

Note added during revision 

During the revision of this paper, [8^ published a preprint in which they pointed out 
that, in the effective model with a time varying curvature, the redshift can be calculated 
from first principles, as it is now done in Section 3.3 of this paper. This refinement of 
the redshift calculation indeed introduces quantitative changes, actually enhancing the 
differences between the model presented here and a FLRW model (as compared to 
the former version of or paper), but does not affect the conclusions of the paper. A 
few comments are in order to clarify the differences between our analysis and the one 
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proposed in [80j. One must note that instead of performing a full MCMC analysis of 
the model, the authors of |80] fitted their model to a ACDM fiducial model. Although 
it leads to correct models that actually reproduce the data (because ACDM is a very 
good fitting model), it is by no way guaranteed that such a procedure provides all the 
acceptable values for the parameters, nor that it provides the most probable ones, as 
clearly shown by our analysis. 
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